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Abstract 

Spherically symmetric solutions in Brans-Dicke theory of relativity with zero 
coupling constant, uj = 0, are derived in the Schwarzschild line-element. The 
solutions are obtained from a cubic transition equation with one small pa- 
rameter. The exterior space-time of one family of solutions is arbitrarily close 
to the exterior Schwarzschild space-time. These nontopological solitons have 
some similarity with soliton stars, and are proposed as candidates for approxi- 
mate black holes for the use in numerical relativity, in particular for treatment 
of horizon boundary conditions. 
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I. INTRODUCTION 



Numerical relativity is considered the prime vehicle for the prediction of gravitational 
wave forms in the late stages of coalescence in binary systems of black holes and neutron 
stars. This effort is complementary to recent analytical approaches for the prediction of 
waveforms in the early stages of coalescence [2,1]. Prediction of waveforms is an important 
factor in the detection of gravitational radiation by LIGO and VIRGO [6,5]. The greatest 
challenges in numerical relativity involve the simulations of multiple orbits in the presence of 
black holes or horizon formation. One of the open problems is proper numerical treatment 
of horizon boundary conditions. 

In this Letter^ we explore the possibility of using approximate black holes (ABHs) for the 
purpose of numerical relativity by introducing a small parameter, e, in a modified setting of 
the physics of the problem. The ABHs to be considered are nontopological solitons in that 
they are produced by a real scalar field in an asymptotically fiat space-time, such that the 
exterior space-time of ABHs approximates the exterior space-time of black holes (BHs) as 
e — > 0. As a first step in this direction, the existence of such solitons shall be demonstrated 
in the Schwarzschild line-element, leaving the general case as an open problem. The fact 
that the scalar field is real, and that the stress-energy tensor in Brans-Dicke does not stem 
from an additive contribution to the Lagrangian ^, however, makes these objects distinct 
from existing soliton stars. 

Concentrating on the problem of black hole coalescence, and thereby ingoring the pres- 
ence of matter, we are at liberty to consider Brans-Dicke theory with = as an approx- 
imation of general relativity. As will be clarified further in what follows, the limit of small 
e, instead of large cu, will serve to approximate general relativity in a suitable way. Vacuum 
with = in Brans-Dicke theory is described by 



Gab 



Rah 



-QabR = ^VaVfcc/) = Sab, 




The author thanks Ira Wasserman for emphasizing this point. 
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where Gab is the Einstein tensor, Rab the Ricci tensor, R = the scalar curvature, gab the 
metric and (j) the real Brans-Dicke scalar field. The stress-energy tensor Sab is trace-free, 
whence □(/> = 0. It thus has the property of a unique coupling to the Einstein tensor 
in Gab = fiSab^ namely = 0, 1 [10]. In a sense. Sab is a four-dimensional, scale free 
generalization of surface tension in three-dimensional continuum mechanics. 

Soliton solutions can be found from their transition equation [11]. In seeking soliton 
solutions with approximate exterior Schwarzschild space-times, we shall work with the 
Schwarzschild line-element 



Earlier analytic studies, by Brans and Dicke [7,3,8] and more recently by Campanelli & 
Lousto [4], have employed different line-elements (including isotropic spatial coordinates). 
The 'black hole' solutions of Campanelli & Lousto [4] are non-standard, in that their solutions 
have infinite horizon surface area, and their exterior space-times are 'funnels' characterized 
by a local minimum of of surfaces of constant r. 



We have Sab = ^'^adb4' = ^{K^b^ ~ ^ab)4''^ where F^^^ are the connections and 
/' = Because of □(/> = 0, we have ^ = e^t/^^F^^ and (/>' = — ^e^~, where 

A is an integration constant, which may be taken to be positive. Because F^^ = 



diag(^,-re-^,-rsin2 6'e-^,^e'^-^), we have ^ = e^^^^F^;, = - !^ - In manifestly 





II. TRANSITION EQUATION 



trace-free form 




(3) 



where e"' = 2 



Einstein's equations 
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T) r It) 

ti^ — -^n 



T + ;^)-;^ = 0' + 2)!J, 



Rq 



(4) 



t 1 



R 



-A/A' 



1 / e"' 



To obtain a differential equation for 7, write (f> = e ^ ip. Upon differentiation, we liave 



re 2 



1 / \/ 

' 2 



(A' + !/')(/) + e 2 ?/)'^ wliile substraction of tlie first and tliird of Einstein' 



equations above gives u' -\- \' = [ru' + 2) 



oA+7 



ru' + 2)^e 2 . Combining tlie last two 



identities gives tp' = ^ru'^e^. N 



•low, re ' 



j-. whence e^{l + rf) = [j-)' 



-e^e^^^ru'. 



Thus, we shall work with 



1 + r-f' 



(5) 



1-rX' = e\l - \ru'e^). 

Now, introduce z = log ^, where Tq is a fiducial mass. Then the combinations rf reduce 
d f 

to which we shall also denote by /'. System (5) becomes the autonomous system of 
ordinary differential equations 



i + \e^+^v' = -1, 



1 _ le^+^)u' 



-1 + + 2e^+"', 



(6) 
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in which V"5'a6 = has been used to replace the second of Einstein's equations. Note that 
solutions to (6) are approximately Schwarzschild, wherever \e'^~^"' << 1. 

Next, consider Z 
Let f{z) = e-^W, so that 



2 



'. System (6) gives A' + i^(e-^-l)-^ = OwithZ' + e-^Z = 0. 



[l-Z)ff' + {l-2Z){f-l)f + 4Z'f 
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with Z' + ^ = 0. Because Z is a strictly monotone function (in the presence of nontrivial Sab 



and wherever the solution is defined), we are at liberty to change the independent variable: 
f{z) = g{Z{zj). Note that £f{z) = 
hence Z{1 - Z)g' + (1 - 2Z)(1 - g)g - ^Z'^g^, that is 



^-^g{Z). We shall denote -^g{Z) also by (/', and 



{{Z - Z^)gy - g^{4Z^ -2Z + 1) = 0. (8) 
Writing h = [Z - Z^)g, we get -(|-)' = Y^ll^^ = ^ + (T^F' Evaluation gives 

= 1 _ 4^ + IZ{1 -Z) = \{Z, + Z)(Z2 - Z), (9) 

where is an integration constant. Consequently, (6) gives the reduced equation u'iz) = 
^Z{z). The zeros Z = —Z\ ~ — e and Z = Z2 ~ 1 — 3e as e becomes small. By the 
aforementioned relation Z' + e^Z, ^(^;) satisfies 

= ^^(^ + ^l)(^-^2), (10) 

which we shall refer to as the the transition equation (c/. the analogous, but second order 
equations in [11]). The solutions to the transition equation are sketched in Figure 1. With 
three zeros in its 'potential' ^Z{Z + Zi){Z — Z2), the transition equation possesses two 
transitions for a given e. With the choices e > and e < 0, therefore, there are altogether 
four different possible transitions, that is, four types of solutions. The solutions to the 
transition equation have a transition from either Z = —Z\ ox Z = Z2 to the exterior 
Schwarzschild limit Z = if e > 0, or from the interior Schwarzschild limit Z = to either 
Z = —Zi OT Z = Z2 if e < 0. These solutions can be regarded to be of exterior and interior 
Schwarzschild type, respectively. The first, exterior Schwarzschild type solutions are the 
soliton solutions of interest to us. 



III. NONTOPOLOGICAL SOLITONS 

To integrate the transition equation, write Z'{-^ — ^^^Zi ( ~'~ z+Zi ) ^ ~ "7-^1 -^2, so 
that Zi^z) follows from the transcedental equation 

I 71^1+^2 7, 7^ 



where K is an integration constant and r]{t) = ^^^[Zi + Z2) ~ 1. By the aforementioned 
relations u' = -Z and Z' + e^Z = 0, can be integrated to give 



which, together with (9) and the transcedental equation (11), completely defines the line 
element. 

The solutions can be treated analytically for the special case e = |, in which case the 
transition equation reads Z' = 4Z{Z^ "4)5 -^1 = -^2 = |- With the choice K = 1, the 
transcedental equation can be inverted analytically to obtain the asymptotically fiat soliton 
solution 



(13) 



The solutions of interest are those with Z ~ as 2; — > 00. For small e, Zi + Z2 ~ 1, ~ 
1, Zi ~ e, Z2 ~ 1 and {Z2 — Z)^^ ~ 1. If e > 0, but small, therefore, ~ Ke~^' , and hence 
r. The approximation to the exterior Schwarzschild solution thus follows: 



l-Ke-^ ' 

„\ -I , Z 1 1 

(14) 

e'^ ~ ~ 1 - Ke-^ = 1 - ^, 

1 + — r ' 

after imposing the relationship Kro = 2M, where M the Schwarzschild mass of the hole. 
Because for e < the transition equation shows soliton solutions with Z ^ Z\ ^^Oas^;— >oo, 
z.e., non-Schwarzschild-like exterior space-times, e > shall be considered henceforth. The 
interior space-time (r < 2M) of the solitons is similarly found to be 

eA ^ 1(^)1/^ 

(15) 

\2MI ' 

Note the extremely rapid decay towards zero of both and e'' as r — > 0, while e^~^ ~ 
^(217)"^' Consequently, </)'(r) ~ "2^7(217)""^ as r — > (taking A = 2M), so that (/> ~ 
^(^)~^. Note, therefore, that the densities \/—g(j) and 

2diag(— + 2,-1,-1,-— )— (16) 
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are finite and regular near tlie origin. We sliall interpret tlie solutions as regular, nontopo- 
logical soliton solutions for Z . The coupling to the metric is through Z — > (/> — > Sah- 

The transition represents an effective horizon in that it represents a thin layer of arbi- 
trarily high redshift. Thus, in the classical sense, the horizon is removed by Sah^ and reduced 
to an effective horizon for sufficiently small e. In particular, with e ~ ep = I I2M , where / is 
the Planck length, the distinction between the effective horizon in the presence of Sah and 
a true horizon in the absence of it can only be discussed in the context of quantum gravity 
for black holes. The new solitons with small e can therefore be regarded as approximations 
to Schwarzschild black holes, in that the difference between the exterior space-time of the 
solitons and that of Schwarzschild black holes is within the realm of quantum gravity for 
e ~ ep. We expect that the present arguments carry through in the case of Kerr black holes, 
but have not pursued this point. 

If Sah is indeed physical, the extreme soliton solutions put Hawking's theorem [9] on 
the non-existence of black holes in BD theory with non-constant (f) in somewhat different 
perspective. Hawking's theorem is purely classical and does not take into account (nor 
prohibit) the possibility encountered here of non-black hole, soliton solutions which are 
classically indistinguishable from black holes, their difference being within the realm of 
quantum gravity. Their interior space-times, however, are dramatically distinct from that 
in standard black holes. In this context, note that a Schwarzschild black hole remains a 
solution of (1), by considering (f) =const, so that, in fact, (1) defines two solutions with like 
exterior space-times. The nontopological solitons, therefore, may be considered bifurcations 
from Schwarzschild black holes in space-time with Sah^ where the bifurcation parameter is e. 

IV. CONCLUSIONS 

The approximation of black holes by extreme soliton solutions offers the advantages of 
a global time-like coordinate and a regularization of horizon surfaces due to an everywhere 
finite (but arbitrary large) redshift (and hence allowing for a great degree of freedom in space- 
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time slicing). In light of this, the solitons are proposed as an Ansatz for treating black holes 
in numerical relativity. Perhaps the a priori most immediate aspects to be addressed are 
(i) the nature and handling of the singularity at the origin and (ii) the persistence of the 
soliton solutions during genuine three-dimensional time-evolutions. It would be of interest 
to consider the effect of a finite coupling constant, uj << 1, on the behavior of (j) near the 
origin. The strong interaction of two solitons in a binary system, furthermore, challenges 
their existence, in that radiation of (j) to infinity may threaten to disperse the solitons. It 
should be appreciated that small e serves to reduce radiation of (j) from the region near the 
effective horizon, in view of the associated high redshift. It would take three-dimensional 
experiments to verify these radiation and stability aspects. 
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FIGURE 1. The solutions to the transition equation Z' = \Z{Z + Zi){Z — Z2) for e > 0. 
The curves shown are those for e = |; for general e > the asymptotes Zi ^ Z2. Two 
nontopological soliton solutions are shown, corresponding to the transition la [Z ~ Z2 to 
Z ~ 0) and lb [Z ~ —Zi to Z ~ 0. The two diverging solutions Ila and lib describe finite 
volume solutions, existing for r < Tq only, while the two constant solutions Ilia [Z = Z2) 
and Illb [Z = —Zi) solutions do not approach the exterior Schwarzschild solution Z = 0. 
The soliton solution la is the nontopological soliton of interest. 

FIGURE 2. The nontopological soliton for the case e = 1/100. A transition layer forms 
at r = To with high redshift as e becomes small. The transition layer serves to replace the 
horizon of the correpsonding Schwarzschild black hole. The radius is normalized with respect 
to 2M. The exterior space-time of the soliton approximates the exterior Schwarzschild space- 
time. 
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Solutions to transition equation for epsiion > 
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